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ON BESSEL'S FUNCTIONS OF THE SECOND KIND. 

By Db. Maxime BdcEEB, Cambridge, Mass. 

In the first volume of the Mathematische Annalen* (1869) Hankel has 
given an investigation of some of the properties of Bessel's functions. This 
paper is probably unknown or inaccessible to many who frequently have to 
make practical use of these functions in physics, so that I have thought it 
worth while to reproduce one of the important results contained in itt accom- 
panying it with a few critical remarks. 

Bessel's equation 

is found (cf. Forsyth, J p. 158) to have as a particular solution, the following 
series, which is convergent for all values of x whatever value n may have : 



J4^x) = 



L2j 



'^=,« (— 1)^ 



^^„/>+j> + l).r(i» + l) 

Here F (n) is the Eulerian integral which reduces to (n — 1) ! when w is a posi- 
tive integer, and which is extended to negative values of n in any legitimate 
way, for instance by means of the formula i"(w + 1) = ?i /^(«).§ 

Thus defined, ^(«) is called a Bessel's function of x of the first kind of 
order n. Since the parameter n enters into the differential equation only 
through its square, J_nix) will also be a solution of the equation, and in general 
the complete primitive of Bessel's equation will have the form AJJ^x) + BJ_„(x). 
If however n is a real integer, the relation J„{x) ^ ( — 1)" J^n{x) holds true (For- 
syth, p. 160) ; so that in this case the above sohation of Bessel's equation is 
no longer general. Confining our attention then, as we may do without loss 
of generality, to positive integral values of n we have the problem : To find in 
this case a solution of Bessel's equation, linearly independent of Jn(x). 

If the parameter of the equation were n — An we should have as two 
independent solutions 

<4-An(«) and t^_„fA»(«)- 

* p. 467. DU OyUnderfnnctionen erster und zweiter Art. 

tl hare confined myself to § 1, although the following sections, especially the last two, are 
of the greatest interest to the physicist as well as to the student of the theory of functions. 

t "A Treatise on Differential Equations." I shall frequently refer to this book by merely 
mentioning the author's name. 

§ For the purpose of this note we may confine our attention to real values of n. This restric- 
tion however is in general undesirable and is not made by Hankel. 

ANNAts Math. VI, 4. 
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Any entire linear combination of these functions will accordingly also be 
a solution. Such a combination is 

Adding to the numerator of this fraction the null expression Ji,(a!) — ( — l)"J'_„{a;) 
wiU not affect it, so that we may write our solution 

2Jn 



If, now, we let in approach zero, we obtain 



dn '^^ ' dn 



as a second solution of Bessel's equation when n is a positive integer. This 
function is called a Bessel's function of the second kind of order n. 
Now we obtain at once by differentiation 

dn 
^dir{_-n+p + l)r' 



dn 



= log 



Jni^) + 



l2j 



"p=_« I \\p 



.t, r{p + 1) 



*^) = -.og[|].^_,W + 



-«p = « (_1)P (x 

p=o7X^+T)[2_ 



dn 



In order to differentiate the reciprocals of r{n + ^ + 1) and i"( — n-\-p-\-\) 
with respect to n we will introduce a new function using a notation similar to 
that of Gauss, and write 



<p{x-) = 



_ d log r{x) 



dx 



This function ^ will evidently have the property 

<p{x) + \ = 4,{x + \); 
so that if a; is a positive integer, 



4:{X) 



+ 



+ ...+i + i + l+^(l), 



and ip (1) has approximately the value — 0.5772. In terms of this function we 
shall have 

'M^ = r{x).4,{x), 



X 


■ JM- 


raj1»p=« (_l)p.^(„,+^ + l) 


a; 


'ij 
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and thus we get 

an 

The expression for the derivative of J-ni^) we might of course get in the 

same way, but this would introduce the factor p, — . "V ■.\ which has for 

positive integral values of n and p the indeterminate form g as long as 
p^n — 1. We will therefore divide our series into two parts at this point, 
getting 



dJ-Jx) _ 
dn 



■log 



■J-n{x) + 



+ 






4r(-n+iP + l)]-' 
dn 






In order to perform the differentiation in the cases which still remain resort 
must be had to a simple formula in the theory of the P function ; namely. 

Fix) . r{l — x) .sirxxK = K ; 
from which we get 

= - sin ( — n -\- p -{- l)7t . r{n — p). 



n-n+p + l) 
and therefore, 



d r 1 

dnyr{—n+p + \) 



= — cos ( — n-\-p + 1)7: . r^{n — p) 

sin( — n+^ + iy . r(n — p) . (p{n — p). 



The value of this expression when n and jp are positive integers and ^ < n is 

_ (_ i)-n+j' + i. r{n —p) = — {— If-P + K Fin —p). 

Thus we have finally performed the differentiation indicated above. The series 
we get may however be written in a more convenient form if in the second (and 
infinite) part the new variable of summation q^^p — n be introduced. This 
gives us 

dJ-n{x) _ __^^Jx-\ J fa) ( ^y.^.[xY'-r' nn-p) \x 
^^^_-log|gJ.^_„(a;) ( 1) ^2j ^to Ai?+1)L2. 



— (— 1)"+' 






L^ 



" "y (-i)«^'(g + i) 

4„ F{n+q + \). F{q + 1) 



'x^^ 
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We thus get finally as the development of our " Bessel's function of the second 
kind of order n " (where n of course represents a positive integer) : 



KJx) = log 



sl-«4(«)— i 



L2j 



'x\ -''i'=»-i r{n —p) fx] 2P 

2j 4o r(-p + i)[2\ 



i 



2] ,t r{n+p + l).r{p + l) [2] 



23) 



I have now completed the account (somewhat fuller, and I hope on that 
account more easy to read, than the original) which I wished to give of § 1 
of the paper of Hankel quoted above.* 

In connection with this subject the following points may be mentioned : — 

(1) That the series we have obtained for '-^^^ and ^!^?i converge un- 
conditionally may be seen at once by forming the quotient of the (p -\- 1) st to 
the/)th term, as this quotient will be seen at once to become and remain as 
small as we please for any finite value of x if we take the number of terms p 
large enough. Nevertheless, we took one step above which still requires justi- 
fication. We assumed, namely, that the derivatives with respect to n of the 
series 



— 2' 



(-IF 



^ZoJ'in+p+l).r{p + l) 



^"•^ ,i'o Ti~n+p+l)J\p+l) 



2 



could be found by differentiating term by term. This, however, requires proof 
as there are many series (for instance, many trigonometric series) which cannot 
be so differentiated. 

Now it is well known that a series may be differentiated term by term if 
the series obtained by performing the differentiation in this way converges 
uniformly for such values of the variable as we areconsidering. We obtained 
from the first of the above-written series by differentiating term by term the 
series : 



p=r (— i)Ps^(w +jg + 1) r«i '^. 



.2J 



and the second of the above series gave us, if we neglect a finite number of 
terms which need not concern us here , 

^ "■> ,L r(n + p + \) r{p + 1)\2\ • 

» The closing part of this section (quoted by Forsyth, p. 165, Ex. 1) of which the above may 
be considered a special case, would seem to be of less general interest. 
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These are the series which must be proved uniformly convergent for all 
positive values of n. This we can do at once by means of the theorem :* 

When the series of the numerically greatest values assumed by the terms 
of the series in the interval in question converges, the series is uniformly con- 
vergent. 

The terms will, however, all become numerically greatest for w ^ 0, and 
for this, as well as for every other value of n, the series have been seen to con- 
verge. 

(2) The function JJx) satisfies for all values of n the fundamental rela- 
tions, which might be called with Gauss relationes inter contiguas : 

nJ„{x) = ^x \Jn+i{x) + J„-i{x) ], 

In fact this seems to be the reasont why the coeflScient of «" in the develop- 
ment of Jn{x) is taken as \2''r{ii -(- 1)]~'; for if it were given a simpler value, 
these relations would not be satisfied. 

It is accordingly of the highest practical importance, that the new func- 
tion KJx) should also satisfy similar relations ; and from our method of intro- 
ducing this function it is evident that it will do so. For we have 

{n — An) Jn-An{x) = i « [^n + l- An(«) + '^n-l-iinix)], 
— {n — Jn) J-n+^nix) = \X [«/_„ + i + dn(«) + J-n-\^■ c.n{«)'\- 

From these equations we get 

^ 2Jh J ■ 

Or, in the limit, when Jn = 0, 

nA'„(a;) = i «![£:„ + i{x) + ^„_i(a;)]. 

The other fundamental relation can be proved in exactly the same way. 

*cf. Harnaek: Die Elemente der Differential- und In tegralreehnung. Leipzig, 1881. Trans- 
lated by G. L. Cathcart, 1891. Book III, Chap. IV, should be consulted. 

1 1 mean the reason which is of weight with practical physicists. Mathematically there are 
several other reasons for this choice, which, however, have nothing to do with the definition of 
Jn(^x) as a solution of a differential equation. 
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(3) We may, however, find another second solution of Bessel's equation sat- 
isfying these relations by adding to KJx) defined as above AJJ(X),viheive A is a 
function neither of n nor of «.* For practical purposes it will be advisable 
to simplify KJ^x) as much as possible by this means. The first term of Knii"^) i^ 
log (aj/2) . Jn{x) which may be reduced in this way to log x . Jn(x). We have 
already noticed that when a; is a positive integer, 

i^(«') = ^ + ^ + . . . + J + i + 1 + <P{1). 

The constant ^(1) may, however, be omitted from all our ^'s ; since it is readily 
seen that by doing so we merely add ^(1) . «/„(«) to Kn{x). We can then write 
our new solution, which I will denote by {K^ix)}, in the following form, hardly 
less adapted to numerical calculation than the ordinary development of JJx) : 



{x^ -"■p=\r\n—p — \)\ 



{K^{x)}= log X JM - i U r" ^" ^^^-fr 

|_^J p=0 jf- 



2p 



-i 



lj":i:(.^i^-^*-^'---l-^-'---^][ir 



(4) Finally it remains to mention the other forms in which the solution 
jr„(a;) is usually found,t namely, a series proceeding either wholly or in part 
according to Bessel's functions of the first kind Jp^x), the positive integer p 
being the variable of summation. On the contrary Hankel's series is merely a 
power series^ except for the necessary logarithmic terms. Hankel's expression 
for A''„(a;) is therefore the most elementary in form, and I think no one who is 
familiar with the methods of obtaining the other expressions for IC„{x) will deny 
that Hankel's expression is also far more easily found. Moreover, his has the 
advantage, as all students of the modern theory of linear differential equations 
know, of being a general method, applicable to the large class of cases where 
logarithmic terms occur in the solution. Why should the working physicist 
not have the advantage of this simple and in the highest degree practical 
theory ? 

• Of course we might also multiply K„(x) by any factor not involving n or x. 

t G. Neumann : Theorie der Besselschen Functionen. 1867. 

Lommel: Studien iiber die Besselschen Functionen. 1868. Also Forsyth, § 105. 

{The function J'„(«) which occurs in the above expression for K„{x) could of course be ex- 
pressed as a power series too, if desired. The above form, however, will probably be as convenient 
in most cases, as K,{x) is rarely used except when Jn(x) has also to be computed. 



